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ABSTRACT 

Recent hydrodynamical simulations have shown that differentially rotating neutron 
stars formed in core-collapse supernovae may develop global non-axisymmetric insta- 
bilities even when T/|iy| (the ratio of the rotational kinetic energy T to the gravita- 
tional potential energy \W\) is relatively small (less than 0.1). Such low-T/|iy| insta- 
bility can give rise to efficient gravitational wave emission from the proto-neutron star. 
We investigate how this instability is affected by magnetic fields using a cylindrical 
stellar model. Wave absorption at the corotation resonance plays an important role 
in facilitating the hydrodynamic low-T/|iy| instability. In the presence of a toroidal 
magnetic field, the corotation resonance is split into two magnetic resonances where 
wave absorptions take place. We show that the toroidal magnetic field suppresses the 
low-T/|VF| instability when the total magnetic energy Wb is of order 0.2 T or larger, 
corresponding to toroidal fields of a few x 10^^ G or stronger. Although poloidal mag- 
netic fields do not influence the instability directly, they can affect the instability 
by generating toroidal fields through linear winding of the initial poloidal field and 
magneto-rotational instability. We show that an initial poloidal field with strength as 
small as 10^"* G may suppress the low-T/jM^I instability. 

Key words: gravitational waves - hydrodynamics - instabilities - MHD - stars: 
neutron - stars: rotation 



1 INTRODUCTION 

Rotating neutron stars (NSs) formed in the core collapse of a massive star or the accretion-induced collapse of a white dwarf 
maybe subject to nonaxisymmetric instabilities (e.g, Andersson 2003; Stergioulas 2003; Ott 2009). The onset and development 
of these rotational instabilities are often parameterized by the ratio P = T/\W\, where T is the rotational kinetic energy and 
W the gravitational potential energy of the star. In particular, the dynamical bar-mode (m = 2) instability sets in when 
P > 0.27 and grows on the dynamical time-scale. This critical /3, originally derived for incompressible Maclaurin spheroid 
in Newtonian gravity (Chandrasekhar 1969), is relatively insensitive to the stiffness of the equation of state as long as the 
degree of differential rotation is not too large (e.g.. Toman et al. 1998; New, Centrella & Tohline 2000; Liu & Lindblom 
2001), although simulations show that it tends to be reduced by general relativity effect as the compactness of the star M/R 
increases (Shibata, Baumgarte & Shapiro 2000; Saijo et al. 2001). There also exist secular instabilities, which are driven by 
some dissipative mechanisms, such as viscosity and gravitational radiation. In the latter case it is known as the Chandrasekhar- 
Friedman-Schutz instability (Chandrasekhar 1970; Friedman & Schutz 1978). Although the threshold of the secular bar-mode 
instability (/3 > 0.14) is easier to be satisfied than the dynamical bar-mode instability, it grows on a much longer time-scale 
due to its dissipative nature (e.g., Lai & Shapiro 1995; Andersson 2003). 

The nonlinear development of the dynamical bar-mode instability has been extensively studied in a large number of 
numerical simulations (Tohline, Durisen & McCoUough 1985; Pickett, Durisen & Davis 1996; Gazes & Tohline 2000; Brown 
2000; Liu 2002; Shibata & Sekiguchi 2005; Camarda et al. 2009). In the early 2000, it was found that for stars with sufficiently 
large differential rotation, dynamical instability can develop at significantly lower /3 than 0.27 (Centrella et al. 2001), even 



* Email: wenfu@astro.cornell.edu (WF); dong@astro.cornell.edu (DL) 



2 W. Fu& D. Lai 



for P on the order of 0.01 (Shibata, Karino & Eriguchi 2002, 2003; Ott et al. 2005; Ou & Tohline 2006; Saijo & Yoshida 
2006; Cedra-Duran, Quilis & Font 2007; Corvino et al. 2010). These low-r/|W| appear to have quite different physical origin 
from the canonical bar- mode instability (see below). Most importantly, recent 3D simulations of a large sample of rotational 
core-collapse models carried out by Dimmelmeier et al. (2008), which include a state-of-the-art treatment of the microphysics 
during collapse and the initial rotational profiles obtained from models of precollapse evolution of massive stellar cores, have 
shown that in many of the models, the proto-NSs exhibit sufficient differential rotation to be subject to the low-T/|W| 
instability (see also Ott et al. 2007) 0. Such proto-NSs would generate strong gravitational waves (GWs), much stronger 
than a non-rotating core-collapse would produce, significantly increasing the possibility of detecting GWs from extra-galactic 
core-collapse supernovae by LIGO and other ground-based GW detectors (Ott 2009) . We note that our current understanding 
of the angular momentum evolution of pre-supernova stars is uncertain, so one cannot predict the rotation profile of the 
collapsing core with great confidence (Heger, Woosley & Spruit 2005). Therefore the detection (or non-detection) of the 
rotational signature of proto-NSs by GW detectors (such as Advanced LIGO) may provide valuable information on massive 
star evolution and the mechanism of core-collapse supernova explosion. 

Despite clear numerical evidence for their existence, the physical origin of the low -r/|Vl^| instabilities remains unclear. 
It has been suggested (Watts, Andersson & Jones 2005; Saijo & Yoshida 2006) that the instabilities are associated with the 
existence of corotation resonance (where the wave pattern speed equals the background fluid rotation rate) inside the star and 
are thus likely to be a subclass of shear instabilities. Corotation resonance has long been known to be the key ingredient for 
some instabilities in other astrophysical fluid systems, such as the Papaloizou-Pringle instability for accretion torii (Papaloizou 
& Pringle 1984; Fu & Lai 2010b) and the corotational instability for thin accretion discs (Narayan, Goldreich & Goodman 
1987; Tsang & Lai 2008, 2009; Lai & Tsang 2009; Fu & Lai 2010a). In addition, numerical calculations by Ou & Tohline 
(2006) suggested that the presence of a local minimum in the radial vortensity profile of the star is also needed to amplify 
the mode (see also Corvino et al. 2010). 

An important issue concerning the low-r/|VK| instability is the effects of magnetic fields. Proto-NSs are expected to 
contain appreciable magnetic fields. In particular, large toroidal fields can be generated from twisting relatively weak poloidal 
fields by differential rotation or from magneto-rotatioinal instabilities (e.g., Balbus & Hawley 1998; Akiyama et al. 2003; 
Obergaulinger et al. 2009). While magnetic fields have a negligible effect on the high r/|H^| instability (Camarda et al. 2009), 
it is not clear whether low-r/|W| can survive in the presence of B fields. Indeed, our previous work on magnetized discs 
showed that even a weak magnetic field can change the structure of corotation resonance significantly (Fu & Lai 2010a). In 
this paper, as a first step of clarifying this issue, we carry out eigenvalue calculation of the effects of purely toroidal B fields 
on low r/|iy| instability by employing a cylindrical stellar model. This paper is the third in our series devoted to study the 
effects of magnetic fields on the global instabilities of various astrophysical fiows, with the previous two focusing on black-hole 
accretion discs (Fu & Lai 2010a) and accretion toil (Fu & Lai 2010b), respectively. 

Our paper is organized as follows. In section 2, we describe the equilibrium model for our rotating magnetized star. In 
section 3, the linearized perturbation equations are presented and boundary conditions derived. In section 4 we present results 
from our numerical calculations. Final summary and discussion of our results are given in section 6. 



2 EQUILIBRIUM MODEL OF A MAGNETIZED ROTATING CYLINDER 

We consider a rotating star with purely toroidal magnetic fields and assume a polytropic equation of state 

P = Kp^ = Kp'+^"', (1) 

where P and p are the gas pressure and density, K, V and N are constants. Although hydromagnetic stellar equilibrium 
models can be constructed easily using the HSCF method (Hachisu 1986; Tomimura & Eriguchi 2005; see Lander & Jones 
2009 for recent works on uniformly rotating stars), linear eigenvalue analysis for such models is difficult. Thus we follow the 
setup in Saijo & Yoshida (2006) by treating the star as an infinite cylinder. We adopt the cylindrical coordinates (r, (j), z). 
All the background variables are assumed to be functions of cylindrical radius r only. The equilibrium state of the cylinder is 
determined by force balance equation in the radial direction 

ldP^ _d£^^^2 _ lrf^__B^^ 
p dr dr p dr Anpr ' 

where Q is the flow rotation rate, is the toroidal magnetic field strength, Pm = B'^/Stv is the magnetic pressure, and "I> is 
the Newtonian gravitational potential which relates to density p via Poisson's equation 

= 4TvGp. (3) 

^ By contrast, the threshold for the canonical bar-mode instability is never reached even when the precollapse core has a very large 
angular momentum, because in that case core bounce would occur at low densities. 
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Eliminating $ from Eqs. ^ and Q yields 

dr (p dr ) 47rGpr + ^^{r Q ) dj. (^p dj- ^ dr ^47rp^ ' 

For numerical convenience, we nondimensionalize variables as follows 
p = ppc = O^pc (5) 



r — r 



{N + l)KpV''-' 

4nG ' ^ ' 



Q = Cl^AnGpc, (7) 



B^^B4,sj4TT{N + l)Kpl+'^'^, (8) 

where pc is the central density and the hatted variables denote dimensionless quantities. Similar dimensionless variables for 
other quantities can be constructed from the list above. We follow Saijo & Yoshida (2006) to adopt the following rotation 
profile: 

" = ^' 

where A and C are constants. The toroidal B field profile we employ is 

B^ = hf{R~f) (10) 

where R denotes the dimensionless boundary of the cylinder and the constant b specifies the field strength. For simplicity, 
we will omit the hats on all variables hereafter unless otherwise noted. The above profile implies that = at both the 
center and the surface of the star. For small r, we have ~ brR, implying a constant axial current for r — > 0. Eq. (|4|) in 
dimensionless form now reads 

d'e , \l ,._iv_i.2 w„ " ^+9'' +Ub\R~r){R~2r)^b'r{2,R-Ar)\e-'' = 2^^{ra). (11) 

dr dr 



dr'^ ^ 



i - Ne-^~^b''r{R -r){R- 2r) 



For an nonmagnetized star, Eq. pi|l reduces to Eq. (3.2) in Saijo & Yoshida (2006). In the limit of zero rotation, it recovers 
the well-known Lane-Emden equation in cylindrical geometry. 

In the limit of B^ = 0, we can simply integrate Eq. (jlip starting from r = using boundary condition that 6=1 and 
^' = to a point where 9 goes to zero, which defines the cylinder surface R. The hydrodynamic equilibrium can thus be easily 
constructed. When is non-zero, we choose an initial guess for the surface radius R based on results for the equivalent 
hydrodynamic model and integrate Eq. imposing the same boundary condition at the center. We stop the integration 
at r = i? to check the value of 9. We then adjust our guess for R and go through the same process, until 9\r^R comes close 
enough to 0. For a given equilibrium state, the rotational kinetic energy T, gravitational potential energy W and magnetic 
energy Wb of the cylinder have the following form 

f pr^Q.'^dV = 



W = - I Pr'^dy = ( / S^rdr 



Bi R%' 



where all the variables are dimensionless and the corresponding physical unit for energy is {N + ifK^p^^l^. Examples of 
the equilibrium density profile are given in Fig. [1] We see that the density profile is not always monotonic: for large C and 
small A (i.e, large rotation rate and large degree of differential rotation), the density maximum is off-centered. 



3 LINEAR PERTURBATION ANALYSIS 
3.1 Perturbation equations 

The cylindrical flow we are considering satisfies the usual ideal MHD equations 

1^ + V ■ {pv) = 0, (13) 
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Figure 1. Density profiles for hydrodynamic and hydromagnetic equilibria. The polytropic index is A'^ = 1. The solid, short-dashed 
and long-dashed lines represent different rotation profiles for non-magnetic models, while the dotted line is for a magnetic model with 
V1/b/|W^| = 0.03. 



^ + V)i; = -ivn-V$-H-^(B-V)B, (14) 
ot p iivp 

^=Vx(^xB), (15) 

V-B = 0, (16) 

= AttGp, (17) 

where 11 — P + Pm is the total pressure. We apply linear perturbations to the above equations by assuming the perturbation of 
any physical variable / to have the form 5f oc e""'*~'"* with m being the azimuthal mode number and lj the wave frequency. 
The resulting linearized perturbation equations contain variables Sv, 5p, 511, 5$ and 5B. To simply the algebra, we define a 
new variable 

Sn 5P B SB 
5h^ — = — + — . 

p p iTTp 

Using Av = Sv + ^ ■ Vv — d^/dt = —ito^ + {v ■ V)^ — —iu^ + nd^/d(j}, we find that the Eulerian perturbation Sv is related 
to the Lagrangian displacement vector ^ hy Sv — — iiD^ — rJl'^i-t^ (prime denotes radial derivative). In terms of ^r, Sh and 5$, 
the MHD perturbation equations (in dimensionless form) can be cast into four first-order differential equations: 

§^=Ai^^r + Ai2Sh + Av,S<i> + Aii^, (18) 
dr dr 

dSh . , . c, A . dS<b 

-7-=A21^r + A22Sh + A2sS^ + A2i—;-, (19) 

dr dr 
dS'^ , ^ , „ A , dS^ 

—;;—= AziS^r + Az2Sh + A-izS^ + A-n—— , (20) 
dr dr 

d f dS^\ . ^ . , (-^ . dS^ 

dr \ IF I " ^^^'^ ^ ^ ^435$ + ^44-^, (21) 
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(22) 



where 

~ 4 2 

(ci +u2_^)((i2 -m2aj2^)((i2 -a;f) r2((I;2 _ ^2^^^) ' 

Ai3 = (24) 

Ai4 = 0, (25) 

4 -2 2 2 4(mcjA</, , d, 2 o2\ , ^ 2 ^2.r dp g dp 

(jj^ — m^Lo\^ dr N V ' p dr p dr 

+ 71— — 2~Tr^ ^ 2 2 \f~-2 2T I'" [('^^'?> " f^^)'^^ + '^i^'^^] + - > (26) 

^ _ rc;}^ [(luj^ - + gji^g;^] gu^ ^ 2m(ma;i^ + cOfl) 1 dp 

(c2+wi^)(i:)2_^2^2^)(j;,2_j^2) (^2 _^ „2 ^)(^2 _ ^2) r((i2 _ rn2a;2 _^) p dr ' 

2mn , , 

A23 = — ^, 28 
rcj 

A24 = -1, (29) 

^31 = A32 = v433 = 0, (30) 

^34 = 1 (31) 

A4. = ~^-^+pf^-lVn-P^-^, (32) 
r dr \ uj^ J ru}^ ru 

A,2=p(!lf4±-l]A.2+p4„ (33) 

4 2 2 

^43 = P^4^ + — , 34 

A44 = -i. (35) 
r 

In the above expressions, lj = uj — mfi is the wave frequency in the co-rotating frame, p = 9^ is the dimensionless density, 
Cs ~ \fdPjdp = \/6/N is the dimensionless sound speed, 

/ 52 

^^<^ = \/ = - r)0-^''' (36) 
y 47rp 

is the toroidal Alfven velocity, ujAtf, ~ VA4>/r — b{R — r)9~^^'^ is the toroidal Alfven frequency. 



c 



2 ,,,2 ■^A4, 



ci + V 



m^<. (37) 



A<t> 

is the slow magnetosonic wave frequency for k — {m/r)(f>, and 

g^^^rn^.^- [h\{R - r){2R - 3r)] (38) 
is the gravitational acceleration in radial direction. 



3.2 Boundary conditions 

To solve Eqs. (|18p -H2i p as an eigenvalue problem, we need four boundary conditions. The outer boundary conditions are 
straightforward. From the perturbed Poisson equation, we know the perturbed potential outside the star scales as 5$ oc r"™. 
By requiring this potential to match smoothly with the potential inside, we obtain our first outer boundary condition: 

^ + II1S$^0 at r = R. (39) 
dr r 

Requiring the Lagrangian pressure perturbation to vanish at the stellar surface yields 
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5h+—ir^Q at r^R. (40) 
dr 

The inner boundary conditions are more involved. As r — > 0, we observe that p — > 0, $7 — > constant, p — >■ constant, ujAtf, oc 
B^/r — >■ constant and 5p is finite. Thus, near the center of the star, Eqs. (|18p -(|2i p can be simplified as 

dCr _ X + mY£_r 5h (5<1> 

dSh X'^ -Y'\ mY Sh 2mQ.S$ dS^ 

-r = — Y — + — + —~ 42 

dr X X r LJ r dr 

1 d /2dS^\ _ Q /^gN 

dr \ dr J r^ ' 
where, 

X = ij'^ — m^ijj\^, Y = 2{mu)\^ + ujil). (44) 
Since perturbation equations are linear, we simply take the solution of Eq. (|43|l to be 

(S$ = r'", at r~0. (45) 
The perturbations and generally take the form 

Cr- = Cir'"-i+C2r'""Mnr, <5/i = Cgr" + C4r'" In r, (46) 

where Ci, C2, C3, C4 are constants so that perturbations remain regular at the center. Eq. (|46|) represents the leading terms 
of the Frobenius expansions of these functions. 

Substituting solution (|46[l into Eqs. (|4ip and (1421) leads to two equations which have the structure ai +02 Inr = 0, where 
Hi, 02 are constants that depend on the values of m, Q., uj, (jja4,, Ci, C2, C3, C4 near the center. Since these equations 
should be satisfied everywhere around the center, we demand ai = a2 = 0. This yields 

C. + ^^^^C.^I^C3 + ^, (47) 

A A UJ'' 

mC4 = {X + Y)C2, (48) 

^ m(x~y)^ x2-y2 2n-cD 

C4 + ^-r^ -C3 = Ci+m z , (49) 

J\ LiJ 

from which we can determine three constants 

m^{m + 2) 2 

<-2 = ^ cj^0, (50) 



,2 



m 



(m + 2)X + -K 2 



C4 = 2 j^^A</>, (51) 



m?{m + 2) 2 m X 

x + r"" ' 2{x + Y) "^^"^^ ^^Ty' 

once we specify C3. When we solve the eigenvalue problem, C3 will be determined together with the eigenfrequency ui. We see 
that since for the specific profile we are considering (see Eq. [10]), ajA0 remains approximately constant near the center, 
C2 and C4 are both finite. Therefore the logarithmic parts in solution (146 |l cannot be neglected. 

In the hydrodynamic limit (B^ = 0, u}a4> = 0), we have C2 = C4 = so that the solutions of and 5h take a purely 
power-law form, and Eq. (|52p reduces to 

= — ^1— (C3 + I), (53) 
uiiui + 2S2) 

which is equivalent to 
Tin 

rLij(ix) + 111) 

Again, constant C3 = 5h/S^ will be determined as a part of the eigenvalue problem. Clearly, for the magnetic cases where 
d In B^/d In r > 1 at r ~ so that cja^ — >■ as r — > 0, the above hydrodynamic boundary condition is also valid. 



3.3 Cowling approximation 

In the Cowling approximation, we neglect the gravitational potential perturbation 5$. The perturbation equations then 
become 



Low-T/\W\ instabilities in proto-neutron stars 7 




T/|W| 



Figure 2. The m = 2 mode frequency as a function of T/IVF] with and with out Cowling approximation. The upper and bottom panels 
show the real and imaginary parts of the frequency, respectively, with Qc being the rotation frequency at the center. The star has no 
magnetic field and the polytropic index is = 1. 



§^ = Aii^r + Ai2Sh, (55) 
ar 

^ = ^21^, + A225h, (56) 
ar 

with the four coefBcients given by the same equations as before. Similarly, for r — )■ the simplified version of Eqs. (jdip and 
are 

d^r _ X + mY S,r Sh 
dSh X'^-Y^^ mYSh 

The outer boundary condition in this case is again given by Eq. (|40p 

Sh +—^r^0 at r^R. (59) 
dr 

The inner boundary condition can be obtained by substituting the power-law solutions oc r™"^ and 5h oc r'" into Eqs. (|57[) 
and (f58)) . giving 

TTl 

Note that for either m = 2 perturbations or a unmagnetized flow, the above inner B.C. reduces to the same form 
ruj(U! + 2il) 



4 NUMERICAL RESULTS 

For most part of this section, we will employ the standard shooting method (Press et al 1992) to solve the two ODEs, Eqs. (I55|l 
and (|56p . subject to boundary conditions (|59p and (|60|l . We focus on the effects of toroidal magnetic fields on the low-T/|VF| 
instability previously found for purely hydrodynamic stars. 
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Figure 3. The m = 2 mode frequency as a function of Wb/|W^| for stellar models with different rotation profiles (thus different T/|Ty ['s). 
The upper and bottom panels show the real and imaginary parts of the frequency, respectively, with Qc being the rotation frequency at 
the center. The other parameters are the same as in Fig. [2] 



0.008 



0.006 



3 0.004 - 



0.003 



— 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



Wb/|W| = 
Wg/lWlc^O.OOl 




0.03 



0.04 0.06 

T/|W| 



0.08 



0.1 



Figure 4. The m = 2 mode growth rate as a function of r/|Vy| with and without toroidal magnetic field. The other parameters are 
A = 0.86 and Af = 1. 
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Figure 5. Example wavefunctions of unstable low-T/|VF| (~ 0.095) mode with A = 0.86, C = 0.9, m = 2 and N = 1. The left column 
shows the radial displacement as a function of radius whereas the right column shows the radial velocity perturbation, with the solid 
and short-dashed lines representing the real and imaginary parts, respectively. The upper and lower panels are for nonmagnetic and 
magnetic stellar models, respectively. The dotted lines indicate the location of the corotation resonance (in the nonmagnetic case) or 
slow magnetosonic resonances (in the magnetic case). The vertical scales of the wavefunctions are arbitrary. 




Figure 6. Angular momentum carried by the wave as a function of r. The model parameters are the same as in Fig. [5] The upper and 
lower panels show the nonmagnetic and magnetic models, respectively. The locations of the corotation resonance and slow resonances 
are indicated by the vertical dotted lines. 
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Before moving on to our main results, let us first examine the validity of Cowling approximation for low-r/|H^| instability. 
To this end, we compare the eigenfrequency calculation with and without Cowling approximation. In Fig. [2] wc fix one of the 
rotation parameters. A, while change the other, C, to obtain different values of r/|VK|. We see that for the whole range of 
T/\W\ considered, the real part of the mode frequency does not show much difference between those two cases. The bottom 
panel of Fig. [2] shows that, for the low r/|W| range (< 0.1), the mode growth rate exhibits qualitatively similar behavior with 
an approximate factor of 2 difference between the two cases. For relatively large r/|ty| (> 0.2), the growth rate no longer 
follows the similar trend when r/|W^| increases. Shibata et al. (2002) also found from their hydrodynamic simulation of a 
similar stellar model that the mode growth rate declines beyond certain T/jVyj. The solid lines ("no Cowling") in our Fig. [2] 
agree well with the results depicted in Fig. 4 of Shibata et al.. With Cowling approximation, we find that the growth always 
increases with increasing r/|W|. Overall, Fig. [2] shows that using the Cowling approximation captures the essential feature 
of the low-T/jM^I instability, especially when T/|VI^| is not much larger than the threshold. 

Figs. [3] and U contain the most important results of this paper. In Fig. [31 we plot the eigenfrequency of the m = 2 
mode as a function of Wb/|VK| (the ratio of magnetic energy to gravitational energy) for different rotation profiles. Note 
that as we change the magnetic field strength, the equilibrium structure, therefore r/|VK| will also change. However, for the 
range of VKb/|VK| we considered, the modification to the equilibrium structure is so small (see dotted line in Fig. [1] for the 
small modification) that r/|W^| is approximately a constant along the three curves. Fig. [3] demonstrates that the low-r/|ty| 
instability can be suppressed by the toroidal magnetic field. The point where the mode growth is completely suppressed 
corresponds to Wb/IM^I ~ 0.2 T/\W\. Fig. [4] shows the mode growth rate as a function of T/|VK| for stellar models with 
different VKb/|H^|. We see that the finite magnetic field shifts the curve towards larger In particular, the magnetic 

field increases the threshold for the instability from r/|Vy| ~ 0.03 for the nonmagnetic model to r/|H^| ~ 0.035 for the 
VKb/|VK| ~ 0.001 model. This finding can be easily understood: increasing rotation drives the instability, whereas magnetic 
field suppresses the instability. Therefore when a finite B field is included, in order to maintain the instability a larger rotation 
rate is needed to overcome the suppressing effect. 

Fig. [5] depicts two example wavefunctions of the overstable low-r/|W| mode. In the nonmagnetic case, the perturbation 
equations are singular at the corotation radius where Cj — uj — mO, — 0. For low-r/|W| modes, the corotation resonance lies 
inside the star, so both the radial displacement and the gradient of the radial velocity perturbation undergo large variations 
across the corotation resonance (see the upper panels). In the magnetic case, however, the corotation resonance is no longer 
a singularity. Instead, the perturbation equations are singular at two slow magnetosonic resonances where tZi — ±tJs, with 
lOb given by Eq. (|37() Q Therefore, the wavefunctions exhibit sudden changes at these two particular locations (see the lower 
panels). This splitting of corotation resonance into two magnetic slow resonances can also be seen in the angular momentum 
flux. In Fig. [6] we show the angular momentum carried by the wave across the star as a function of radius (see Fu & Lai 2010a 
for the fiux formula). In the upper panel (nonmagnetic case), we see that F{r) experiences a sudden jump at the corotation 
resonance, whereas in the lower panel (magnetic case), two jumps occur at the two slow resonances and have different signs. 
This is similar to thin accretion discs studied in Fu & Lai (2010a). However, since the outer boundary condition we employed 
here (free surface) is totally different from the one used in Fu & Lai (2010a) (outgoing waves), we cannot directly relate the 
flux jump (or the net jump in the case of two resonances) to the magnitude of the growth rate. In any case, it is clear from 
Fig. [6] that the corotation resonance indeed plays an important role in driving the hydrodynamic low-r/|Vy| instability and 
the toroidal magnetic field affects the instability by splitting the corotation resonance into two magnetic slow resonances. 
The property of the unstable mode in the presence of a magnetic field is determined by the combined effects from both slow 
resonances. 



5 DISCUSSION 

Recent studies of rotating (but nonmagnetic) core-collapse supernovae (e.g., Dimmelmeier et al. 2008) have demonstrated 
that newly formed neutron stars can develop nonaxisymmetric global instabilities with low r/|M^|, and such instabilities 
lead to significant gravitational wave emission. In this paper, we have carried out the linear stability analysis of magnetic, 
differentially rotating stars (modeled as a cylinder) to examine how magnegic fields affect the low-r/|VK| rotational instability. 
We show that the wave absorption at the corotation resonance plays an important role in the instability. In the presence of 
a toroidal magnetic field, the corotation resonance is split into two magnetic resonances, where wave absorptions of opposite 
signs take place. Our main result is that toroidal magnetic fields reduce the growth rate of the low-r/|W| instability and 
increase the threshold r/|W| value above which the instability occurs. To significantly affect the instability, the required 

^ Note that ui'^ — m?uP^^ = is not a singularity even though it appears to be a singular term similar to ili^ — uJ^ in those coefficients 
of differential equations. This apparent singular term, one can show, will be canceled by some subtle mathematical manipulations. But 
this cancelation only works for the particular setup we consider here (pure toroidal B field, no vertical structure in perturbations). In 
general (i.e., with mixed B field or finite fcz), equations will be singular at both u}'^ = mPuj'^, and lu^ = ui'^ (see Fu & Lai 2010a). 
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VKb/|VK| (the ratio of the magnetic energy Wb to the gravitational potential energy \W\) should be of order 0.2r/|W^i or 
larger (see Figs. I3ll4)) . As the critical T/\W\ ranges from 0.01 to 0.1, the required Wb/|W^| lies between 0.002 and 0.02. Using 
\W\ ~ {3/5)GM^/R and Wb ~ (-B^/87r)(47ri?V3), we have 



Wb If Y f R Y f M 



W\ 300 V2x 1016 G/ V20km/ \1AM, 



(62) 



Thus, only toroidal magnetic fields stronger than 2 x 10^® G can significantly affect the low-r/|VK| instability. 

In our simple (cylindrical) stellar model, poloidal (vertical) magnetic fields do not directly affect the rotational instability 
because the unstable modes do not have vertical structure (i.e., the vertical wavenumber is zero). We believe that this also 
holds for more realistic stellar models, although more investigations are needed. 

Nevertheless, even a relatively weak poloidal magnetic field present in the proto-neutron star may indirectly affect the 
T/\W\ instability. In the core-collapse supernova scenario, differential rotation naturally arises inside the stellar core during 
the collapse (e.g., Akiyama & Wheeler 2005; Ott et al. 2006). This differential rotation can generate significant toroidal 
magnetic fields by winding the initial poloidal field and by magneto-rotational instability (Akiyama et al. 2003; Obergaulinger 
et al. 2009). Consider first the linear winding of the poloidal field Bp. The toroidal field grows in time as Btf, ~ BpAflt, where 
Af2 is the difference in the rotation rate across the proto-neutron star. Thus the ratio of magnetic energy Wb ~ B^R? /Q 
and the rotational energy T ~ 0.2MR?{1\QY increases as Wb/T ~ B^Rt^/M. The time to reach a given Wb/T = / is then 
ttwist ~ {fM/B^RY^ = YT/SR/vAp, where VAp = Bp/y/A-Kp is the Alfven speed associated with Bp. On the other hand, 
the growth time of the low-r/|W^| instability is tgrow ~ l/i^i ~ {R^ /GMY^'^ where oji is the dimensionless growth rate 
in units of the Keplerian frequency (oji is approximately the vertical axis of the bottom panel in Fig. |3]). For the low-T/j W| 
instability to operate before being suppressed by the large B^ (generated by twisting the initial poloidal field Bp), we require 
^twist ^ ^grow , i.e., 

, (; 2^) . , . o, (X) - (^) (^) (^) - ^ 

Since Cji is of order 10^'^ or larger, and the toroidal field suppresses the instability when / — Wb/T ~ 0.2 (see Fig. [3|, we 
see that an initial poloidal field strong than 10^'* G can lead to the suppression of the instability. In other words, when the 
initial poloidal field is less than 10^* G, the toroidal field will not grow fast enough by linear winding so that the low-T/|VF| 
instability still has a chance to develop. 

The effect of magneto-rotational instability (MRI) is harder to quantify. In the linear regime, MRI operates in modes 
with vertical structure (i.e., finite vertical wave number), which are independent from the T/| Vl^|-unstable modes studied in 
this paper. However, the nonlinear development of MRI may generate significant magnetic fields (both poloidal and toroidal) 
on a short timescale (of order the rotation period). There have been many MHD simulations of core collapse supernovae (e.g., 
Ardeljan et al. 2000, 2005; Kotake et al. 2004; Yamada & Sawai 2004; Obergaulinger et al. 2006; Burrows et al. 2007). Most of 
these simulations cannot resolve the MRI unless they employ drastically strong initial fields. It has been suggested that when 
MRI saturates, a large fraction of the kinetic energy in the differential rotation is converted to the magnetic energy (Akiyama 
et al. 2003; Obergaulinger et al. 2009), i.e., Wb/T may approach unity on a dynamical time. Our result in this paper shows 
that the r/|W| instability is strongly reduced when VKs/lrj reaches 0.2. Therefore it would be important to quantify the 
saturation field of the MRI in proto-neutron stars. In addition, the MRI can lead to efficient angular momentum transport 
in different region of the star. This may also affect the r/|W| instability. Clearly, these issues must be resolved in order to 
evaluate whether the low-r/|VK| instability can develop in astrophysically realistic proto-neutron stars. 
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